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We combine the QSO samples from the 2dF QSO Redshift Survey (2QZ) and 
the 2dF-SDSS LRG and QSO Survey (2SLAQ) in order to investigate the cluster- 
ing of z ~ 1.4 QSOs and measure the correlation function (£). The clustering signal 
in redshift-space and projected along the sky direction is similar to that previously 
obtained from the 2QZ sample alone. By fitting functional forms to £(<7, 7r), the corre- 
lation function measured along and across the line of sight, we find, as expected, that 
/3, the dynamical infall parameter and fl m , the cosmological density parameter, are 
degenerate. However, this degeneracy can be lifted by using linear theory predictions 
under different cosmological scenarios. Using the combination of the 2QZ and 2SLAQ 
QSO data, we obtain: (3qso(z = 1-4) = 0.601^, VL° m = 0.25^ ^ which imply a value 
for the QSO bias, b(z = 1.4) = 1.5 ± 0.2. 

The combination of the 2QZ with the fainter 2SLAQ QSO sample further reveals 
that QSO clustering does not depend strongly on luminosity at fixed redshift. This 
result is inconsistent with the expectation of simple 'high peaks' biasing models where 
more luminous, rare QSOs are assumed to inhabit higher mass haloes. The data are 
more consistent with models which predict that QSOs of different luminosities reside 
in haloes of similar mass. By assuming ellipsoidal models for the collapse of density 
perturbations, we estimate the mass of the dark matter haloes which the QSOs inhabit. 
We find that halo mass does not evolve strongly with redshift nor depend on QSO 
luminosity. Assuming a range of relations which relate halo to black hole mass we 
investigate how black hole mass correlates with luminosity and redshift and ascertain 
the relation between Eddington efficiency and black hole mass. Our results suggest 
that QSOs of different luminosities may contain black holes of similar mass. 

Key words: 
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1 INTRODUCTION 

There is a significant amount of observational evidence for 
the existence of supermassive black holes in the centre of 
galactic haloes. This conclusion is based on studies which 
span a wide redshift-range. Whilst at low-z, the evidence for 
the presence of black holes comes from dynamical surveys 
of galaxies in the local Universe llKormendv fc Richstond 
ll995l ; lRichstone et al.l ll998; Magorrian et al.lll998l ), at high- 
z, black hole - host galaxy studies are pursued by us- 
ing the width of quasar (QSO) broad emission lines to 
estimate black hole masses and the host galaxy's narrow 
emission lines to d etermine stellar velocity dispersion (e.g. 
IShields et alj|2006h . These results hint at a correlation be- 
tween the growth/physics of the bulge and dark matter 
halo, and the physics of accretion of mas s onto the cen- 
tral b lack hole and subsequent growth (e.g. iTremaine et al.l 
2002). The relation between the bulge and its black hole 
is the s ubject of intense observationa l and theoretical in- 
teres t (iKauffmann fc HaehneltJ |2000|; iFerrarese fc Merritt 



ICroom et ail |2004) data have allowed studies such as 



those of 



Porciani et ail (|2004h . IMvers et all l|2006h . and 



2000; 



Gebhardt et al.ll200d ; lFerraresd[2002l ; IWvithe fc Loeb 



2005b). Many uncertainties still exist when trying to inter- 
pret this black hole - bulge connection. One possible sce- 
nario is that the mechanism that "feeds" black hole growth 
is the same, or is correlated to, those properties responsi- 
ble for bulge growth, such as mergers or instabilities, which 
may also lead to enhanced star formation; some of the gas 
may instead "fuel" the black hole, an d consequently lead to 
QSO activity (e.g. iBower et al.l 1200^ ) . This picture is sup- 
ported by the similar "shape" of the cosmological star for- 
mation history of the Universe and the evolut ion of the QSO 



Bovle et al. 


1988; 


Schmidt et al.ll 19951 


Dunlop et al.ll200S 


)• 



In the standard scenario, QSO activity is triggered 
by accretion onto a sup ermassive black hole (SMBH, 
e.g. IHopkins et alj 120061 ) . Given that the growth of 
the SMBH relates to that of the underlying dark 



matter halo dBaes et all 120031 ; IWvithe fc Loebl l2005al ; 
IWvithe fc Padmanabhan 20061 ') and the halo properties are 
correlated with the local density contrast, clustering mea- 
surements provide an insight into QSO and black hole 
physics. 

QSO clustering measurements allow determinations of 
halo masses and how they relate to black hole mass. QSO 
lifetimes, which have been the basis of interpre tations of 
QSO luminosity functions (IHopkins et alj l2005al) can also 
be in ferred from clustering measurements (e.g. Croom et al.l 
120051 ). and hence permit discrimination between QSO evolu- 
tionary m odels, such as a c os mologically long-liv ed popula- 
tion (e.g. iBovle et al.l 120001 1. iMiller et all (|2006l ) addressed 
the change of accretion efficiency with redshift, arguing 
that, even though the mass of the black holes grows with 
time as galaxies grow hierarchically, the mean accretion 
rate decreases with decreasing redshift, hence leading to 
a decrease of the QSO luminosity with time. This pic- 
ture is supported by theoret ical models, such as that of 
IKauffmann fc Haehneltl (^OOOh . 

The evolution of QSO clustering has been the sub- 
ject of recent studies. In particular, the wealth of infor- 
mation containe d in the Sloan Digital Sky Survey (SDSS; 
I York et al.|[2000l ) and the 2dF QSO Redshift Survey (2QZ; 



ICroom et all rt2005h . who measured the redshift depend' 
of QSO clustering. In particular, the latter inferred the evo- 
lution of halo mass with redshift, besides estimating black 
hole masses and accretion efficiencies, based on QSO clus- 
tering measurements from the 2QZ sample. However, and 
as pointed out by those authors, these studies do not take 
into account any potential luminosity dependence of QSO 
clustering. 

It is not trivial to address the possible dependence of 
QSO clustering on luminosity. Due to the evolution of the 
QSO luminosity function and the flux-limited nature of the 
2QZ and most other surveys, the most luminous QSOs lie 
at high redshifts, while the faintest ones have low redshifts. 
The lowest and highest redshift objects in the 2QZ sample 
extend throughout separate luminosity ranges, hence ham- 
pering any attempt to study the effects of luminosity on 
QSO clustering, black hole masses and accretion efficiencies, 
free from any possible evolutionary biases. 

This necessary caveat in any study of luminosity depen- 
dence of QSO clustering was one of the main motivations for 
the 2SLAQ (2dF-SDSS LRG and QSO) QSO survey. Using 
faint, photometric QSO candidates from the SDSS QSO sur- 
vey, the observations at the 2dF facility result in an exten- 
sion of the previous 2QZ survey to fainter magnitudes. The 
faint magnitude limit of g = 21.85 is ~ 1 magnitude fainter 
than that of the 2QZ, and the new data, spanning a similar 
z-range as the 2QZ, constitute a new, potentially powerful 
tool to disentangle the effects of luminosity and redshift on 
the clustering of QSOs, thus providing a new test of cur- 
rent QSO, black hole and bias models. With its fainter QSO 
m agnitude limit than the photometric SDSS QSO catalogue 
of IMvers et al.1 (|2006h . 2SLAQ, despite its smaller statistical 
weight, constitutes a more valuable tool for breaking the L-z 
degeneracy. 

In this paper we combine the 2QZ and 2SLAQ QSO 
samples and analyse the clustering of z ~ 1.5 QSOs. In 
addition, we use the wide luminosity range covered by the 
combination of the two ensembles to determine the lumi- 
nosity dependence of QSO clustering, free from evolution- 
ary effects. In section [2] we present a brief description of the 
2SLAQ QSO survey. We then measure the clustering signal 
of the QSOs, in redshift-space (z-space); projected along the 
sky direction (and hence free of dynamical distortions); and 
in orthogonal directions (section [3J). These measurements 
allow us to model the anisotropies due to dynamical and ge- 
ometrical distortions in the clustering signal and constrain 

and Pqso(z = 1.5). This analysis is discussed in section 
[4] In section [5] we describe the L-z degeneracy and how we 
attempt to break it by combining the 2QZ and 2SLAQ QSO 
samples. Our QSO clustering measurements as a function of 
magnitude and redshift follow in section[6] We then attempt 
to determine if QSO bias correlates with QSO luminosity, 
and how these results affect the average mass of the dark 
matter haloes the QSO inhabit (section [7}. Assuming that 
the mass of the dark matter halo correlates with that of 
the black hole associated to the QSO, we determine how 
the black hole mass changes with redshift and luminosity, 
and discuss how our results affect the black hole accretion 
efficiency, in section [8] Finally, in section [9] we outline the 
conclusions of this paper. 
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Objects Priority 



Guide stars 9 

Main sample LRGs, sparsely sampled 8 

Remaining main sample LRGs 7 

g > 20.5 QSOs, sparsely sampled 6 

Remaining g > 20.5 QSOs 5 

Extra LRGs and high-2 QSOs 4 

g < 20.5 QSOs 3 

Previously observed objects with good id 1 



Table 1. 2dF priorities. Objects with higher priorities have a 
higher likelihood of being assigned a 2dF fibre. 



2 THE 2SLAQ QSO SURVEY 

The 2SLAQ QSO survey is an extension of the previous 
2QZ survey to fainter magnitudes. The m ain aspects and 
descri ption of this survey can be found in iRichards et al.l 
(2005), who report on the first 3 semesters of the data collec- 
tion and present luminosity function results from the sam- 
ple of ~ 5600 QSOs obtained at the time. Now that the 
survey has been completed and the analysis of the data is 
being developed, there are a total of ~ 9000 (z< 3) QSOs. 
Both the imaging and spectroscopic data, obtained from the 
Sloan telescope and the AAT respectively, are extensively 
described by Croom et al., (in prep.). 

The sky regions surveyed by the 2dF instrument consist 
of two 2° - wide equatorial strips, containing the QSO can- 
didates observed by SDSS survey. Not all of the full strips 
were observed, but rather "sections" of them. Fig. [1] shows 
the two strips, on the NGC and SGC. The NGC photomet- 
ric candidates are shown in green and the SGC ones in pink. 
The blue (red) circles are all the spectroscopically identified 
QSOs in the NGC (SGC). The 2dF pointings are shown as 
black circles. The "sections" in the NGC were indexed "a, b, 
c, d, e" and the one in the SGC "s". Each 2dF pointing was 
labelled with the index of the region where it fell followed 
by a number, which refers to its position along the strip. 

The QSO observations were performed simultaneously 
with those of the LRGs. 200 2dF fibres wer e allocated to 
the L RGs and 200 to the QSO observations l|Cannon et al.l 
l200fj ). The LRG fibres then link to the 2dF "red spectro- 
graph" and the QSO fibres to the "blue spectrograph". Each 
block of 10 consecutive fibres along the edge of the 2dF field 
connects to a different spectrograph, alternately blue and 
red. Therefore, the QSO completeness mask in each 2dF 
pointing shows a "dented structure" along the edge of the 
field, due to the fact that the fibres a re limited to an angle 
of 14° (see, e.g. IRichards et~al1l2005l ). The probability of a 
given QSO/LRG candidate being assigned a 2dF fibre de- 
pends on its priority. The assigned priorities of the objects 
in the input catalogue (see table [TJ will affect the likelihood 
that those objects will be observed. Objects with higher pri- 
ority will have a higher likelihood to be assigned a 2dF fibre. 

Tables [2] and |3 show the number of QSOs, narrow emis- 
sion line galaxies (NELGs) and stars that were observed. 
Ql and Q2 refer to the identification quality: Ql are objects 
with good identification quality and Q2 refer t o objects with 
lower identification quality (see section 2.3 of ICroom et al.l 
120041 for further details on quality identification flags). Over- 



ID 


All 


Ql 


Q2 


QSOs 


6680 (57.89%) 


6482 (56.17%) 


198 (1.72%) 


NELGs 


2077 (18.00%) 


2043 (17.71%) 


34 (0.29%) 


stars 


1829 (15.85%) 


1604 (13.90%) 


225 (1.95%) 


TOTAL 


10586 (92.20%) 


10129 (88.15%) 


457 (4.05%) 


Table 2. Number of QSOs in 


the NGC 2SLAQ strip. 


ID 


All 


Ql 


Q2 


QSOs 


2378(49.68%) 


2282(47.67%) 


96(2.01%) 


NELGs 


905(18.91%) 


881(18.40%) 


24(0.50%) 


stars 


835(17.44%) 


739(15.44%) 


96(2.01%) 


TOTAL 


4118(86.02%) 


3902(81.51%) 


216(4.51%) 



Table 3. Number of QSOs in the SGC 2SLAQ strip. 



all, the sky density of QSO candidates in 138.4 deg~ 2 and 
that of confirmed QSOs is 44.7 deg -2 . 

As we are observing faint QSOs, we also expect them 
to have a higher space density than that achieved from 
other, previous surveys, such as the 2QZ or the SDSS. 
This is evident from the wedge plot in Fig. [2] which 
shows the radial projection of the 2SLAQ NGC strip (in 
pink). The QSOs observed from the 2QZ and SDSS DR4 
ijAdelman-McCarthv et al]|2006l ) datasets are also shown in 
the wedge plot, in blue and cyan squares, respectively. 



3 QSO CLUSTERING 

Completeness issues within a 2dF pointing must be taken 
into account when constructing the angular mask used to 
generate a random set of points, which is necessary to mea- 
sure QSO clustering from the 2SLAQ survey. The complete- 
ness in each pointing depends on two factors: (i) the coverage 
completeness, given by the fraction of QSO candidates that 
were assigned a 2dF fibre; and (ii) the spectroscopic com- 
pleteness, representing the fraction of observed candidates 
which have good redshift quality. In addition, one needs to 
calculate the excess probability of finding a QSO in overlap- 
ping 2dF pointings. 

The fact that the 2dF instrument cannot place two fi- 
bres any closer than ~ 30 arcsec means that an additional 
incompleteness can potentially lead to an artificial deficit of 
close QSO pairs in 2dF surveys. To make an approximate 
correction for these effects, on e can measure the angular 
correlation function, w{9) fe.g. lHawkins et al.ll2003l ). Com- 
paring this to the angular correlation measured in the total 
input catalogue allows one to estimate the average deficit 
of close pair s at s mall angular separations. As shown by 
ICroom et al.l (|200ll ). this deficit is negligible in the 2QZ sam- 
ple. In the 2SLAQ sample, however, the deficit of pairs can, 
potentially, constitute a bigger bias. This is due to the fact 
that, in contrast to what happens in the 2QZ survey, the 
2SLAQ QSOs are assigned a lower observational priority 
than the main sample LRGs. Therefore, the QSO-assigned 
fibres will only be positioned in areas allowed by the under- 
lying angular distribution of the LRG fibres. Fig. [3] shows 
the w(6) measurements of the 2QZ+2SLAQ sample and the 
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Figure 1. A sky map of the 2SLAQ observations. The upper two panels show the NGC, the lower two panels form the SGC. The black 
circles represent the 2dF fields observed. Green and pink dots are the NGC and SGC QSO candidates, respectively. The (blue and red) 
points represent the positions of the (NGC and SGC) spectroscopically confirmed QSOs. 
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Figure 2. This shows a wedge of the sky region covered by the 2SLAQ QSO survey. The other two QSO surveys that have partial 
coverage in the same area are shown for comparison. The number in brackets gives the total number of QSOs with z < 3 per survey in 
this region. The comoving distance is computed assuming a Q 1 ^ = 0.3, = 0.7, h = 0.7 cosmology. 



2SLAQ and 2QZ samples separately. In order to better dis- 
tinguish between the errorbars, the 2SLAQ values are offset 
by a shift of A0 = +0.02 and the 2QZ w(0) points by a 
shift of —0.02. To account for the fibre-collision effects in 
the clustering of the 2SLAQ QSOs, we followed the method 
applied in previous work to the Two-degree Field G alaxy 
Redshift (2dFGRS) survey data ijHawkins et al.ll2003t ): the 



number of QSO pairs at a given separation is assigned a 
weight that depends on the QSO's angular separation. Since 
the QSO sample spans a wide redshift range, the input cat- 
alogue is expected to show zero correlation at all angular 
separations w(0) ~ 0, V0. In this cas e, the weight assigned 
to each QSO pair using the method of lHawkins et al] (|2003h 
is 1/(1 + w(0)). The "imprint" of the LRG angular distribu- 
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Figure 3. The angular correlation function measured for the 2QZ 
survey (solid blue stars and solid blue line), the 2SLAQ QSO sur- 
vey (open blue stars and dashed blue line) and the 2QZ and 
2SLAQ QSO surveys combined (red circles and line). The w(9) 
measurements are very similar in both cases and show that the 
deficit of pairs seen at the smallest scales is not significant at typ- 
ical QSO-QSO comoving separations. Note that the 2QZ values 
are offset by a shift of Ad = —0.02 and the 2SLAQ values by a 
shift of A6 = +0.02. 



tion on the QSO fibres, due to these having been assigned a 
low 2dF priority, is also accounted for: when generating the 
random catalogue for the determination of the correlation 
functions of the 2SLAQ QSO sample, any random point has 
a zero probability of lying closer than 30 arcsec to any ob- 
served LRG. Although these effects have been considered, 
they have negligible effect on our clustering results. 

Equally as relevant is the radial completeness, which 
also needs to be accurately described by the unclustered, or 
"random" distribution. Fig. [3] shows the (0.3 < z < 2.9) 
redshift distribution of the 2QZ and 2SLAQ QSOs, in Az = 
0.13 bins. The red line represents the 2SLAQ NGC while 
the blue line the 2SLAQ SGC. The green and pink lines are 
the z-distributions of the 2QZ NGC and 2QZ SGC QSOs, 
respectively. Dashed lines also show the polynomial fits to 
those distributions that were used to generate the random 
distribution. 

The 2QZ survey comprises 22416 (id quality 1) QSOs 
in the redshift range 0.3 < z < 2.9 (9982 in the NGC and 
12434 in the SGC). The 2SLAQ QSO sample, when impos- 
ing faint magnitude cuts (20.5 < g < 21.85) in addition to 
these z-cuts, comprises a total of 6374 QSOs (4574 in the 
NGC and 1800 in the SGC). The fact that the 2SLAQ N(z) 
is steeper, at low-2, is possibly due to QSO contamination 
by host galaxies, affecting the colour selection of fainter 
QSOs. The median redshift of the 2QZ+2SLAQ sample is 
< z >= 1.50. 

After generating a random catalogue we can then com- 
bine the new 2SLAQ QSO sample with the 2QZ sample, 
and compute the QSO clustering by means of correlation 
functions. We start by estimating £(s), the 2-point correla- 
tion function measured in z-space. This is presented in Fig. 
[5] (f illed red c i rcles) . The estimator used to measure £(s) is 
the lHamiltonl (Il993l ) estimator : 



Figure 4. 2SLAQ QSO and 2QZ N(z). Red line is the NGC and 
the blue line the SGC. The green line represents the 2QZ NGC 
and the pink line the 2QZ SGC. Also shown, as dashed lines, are 
the polynomial fits that were used to model the radial distribution 
of the random points. 



where < DD(s) >, < DR(s) >, < RR(s) > are the mean 
number of QSO-QSO, QSO-random and random-random 
pairs at separation s. For comp arison, also shown is th e pre- 
viously determined 2QZ £(s) |da Angela et all 120051 ). the 
2SLAQ QSO £(s) an d also the £(s) m easurements of the 
2SLAQ LRG sample ijRoss et alJ 120061 ). To make the plot 
clearer, we have offset the 2QZ and 2SLAQ f(s) points by a 
factor of 0.02 and —0.02, respectively. 

Including the 2SLAQ QSO sample does not affect the 
shape of the previously measured 2QZ £(s). The £(s) mea- 
sured from both samples, including or not the 2SLAQ QSOs, 
are indeed very similar. We have verified the statistical 
weight of including the 2SLAQ sample by comparing the 
number of QSO-QSO pairs at separations < 20/i _1 Mpc, 
and verified that the combined 2QZ+2SLAQ sample has 
~ 65% more QSO-QSO pairs within 20/i _1 Mpc than the 
2QZ sample alone. This gain also includes the contribution 
of the cross pairs between the 2SLAQ and 2QZ samples, on 
the NGC strip. The 2SLAQ LRGs have a higher clustering 
amplitude than the 2SLAQ QSOs. At smaller scales the two 
samples also differ in the shape of their correlation functions. 
This is probably due to the different z-space distortions that 
affect the LRGs and the 2QZ and 2SLAQ QSOs, a contribut- 
ing factor to which will be the higher redshift errors of the 
QSOs. 

Also shown are two d ifferen t 2QZ £(s) models, ob- 
tained by Ida Angela et all (|2005l ). The dashed line is the 
best fitting 2QZ power-law model, in the range 5 < s < 
50/!, _1 Mpc (£(s) = (s/6.50)" 1 ' 89 ), and the solid line is the 
£(s) model obtained from convolving a double power-law 
£(r) model (Eq. [5} with the z-space distortions parame- 
terised by < w 2 z > x f 2 = 800 kms" 1 and (3{z) = 0.32. 



f(r) = 



(r/6.00)" 
(r/7.25)" 



r < 10/i _1 Mpc 
r > 10/i _1 Mpc 



(2) 



< DD(s) >< RR(s) > 
< DR(s) > 2 



(1) 



It can be seen that the model is still a good description 
of the joint QSO £(s) measurements, indicating that the 
2SLAQ QSOs should have a similar real-space clustering 
and be subjected to the same dynamical distortions as the 
2QZ QSOs. The fitting of these models does not take into 
consideration the correlations between the errors at different 
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Figure 5. The red circles show the measured from the 

2SLAQ and 2QZ samples and the blue triangles the 2QZ re- 
sults (from da Angela et al., 2005). The £(s) measurements 
are very similar, both in amplitude and shape. The green stars 
show the 2SLAQ LRG measurements (I Ross et al. 2006). The 
dashed and solid lines show two models: the best fitting 2QZ 
5 < s < 50 h" 1 Mpc power law (dashed); and the double power 
law £(r) model, "distorted" by dynamical motions parameterised 
by < w\ >V2 = 800 km s -1 and /3(z) = 0.32. Note that the £(s) 
values from the individual 2SLAQ and 2QZ samples have been 
offset by log As of 0.02 and —0.02, respectively. 



separations. However, Ida Angela! (|2006t ) showed that taking 
into account the full covariance matrix when fitting the 2QZ 
£(s) does not affect the so and 7 values by more than la. 

The errors shown in Fig. [5] are "jacknife" estimates, es- 
timated by splitting the 2QZ+2SLAQ sample in 16 sub- 
samples. We compared the jacknife and Poisson error esti- 
mates in our £(s) computation. The Poisson error estimates 
should, in principle, provide a fair description of t he un- 
certa i nty for the 2QZ QSO c lustering measurements (|Hovlel 
2000; Ida Angela et al.l l2005h . Here we test this hypothesis 
for the new sample containing the 2QZ and 2SLAQ QSOs. 
We divide up the overall 2QZ+2SLAQ dataset into 16 sub- 
samples and compute £(s) in the overall set minus each of 
the 16 subsamples in turrQ. The 16 measurements of £(s) 
are then c ombined as follows , in order to obtain the jacknife 
error (e.g. iMvers et alj|2005h : 



•^"jacknife 



E 



DRi(s) 



•te(s)-6o t (s)) 2 



(3) 



where N is the total number of subsamples (16, in this 

1 This £ (s) computation w as performed using the fcd-tree algo- 
rithm of lMoore et al.l 11200111 . 




Figure 6. Red circles and solid line show the ratio jacknife and 
Poisson £(s) errors. Poisson errors seem to under-predict the un- 
certainty in §(s) at all scales, and considerably at the largest 
scales. At intermediate, 4< s< 20 /i -1 Mpc scales, the ratio of the 
two error estimates is approximately constant and ~ 1.25 (dotted 
line). 



case); the subscript i refers to the whole dataset minus 
subsample i; and tot refers to the whole 2QZ+2SLAQ 
QSO sample. The "DR ratio" accounts for the fact that 
the subsamples may not necessarily contain exactly the 
same number of QSOs. Fig [6] shows the ratio between the 
jacknife and the Poisson errors. It can be seen that, on all 
scales, Poisson errors underestimate the uncertainty on the 
clustering measurements, especially at the largest scales. 
On scales 2< s< 4 /i _1 Mpc, the two estimates are quite 
similar, but on 4< s< 20 /i _1 Mpc scales, where most of the 
clustering signal is obtained, the jacknife errors are, on 
average, 1.25 times bigger than Poisson errors (dotted line). 
At larger scales, where there are fewer QSO independent 
pairs, the Poisson estimates largely under-predict the 
true error estimate a s has been previou sly discussed, (e.g. 
IShanks fc Bovlell 19941 ; IMvers et al.ll2005h . 

Fig. [7] shows the projected correlation function 
measured from the 2QZ+2SLAQ sample (red circles). 
This is very similar t o the previous 2QZ measurement 
(|da Angela et al] 120051 . blue triangles, offset by a factor 
of log Act = 0.02). The open blue triangles represent the 
w p (a)/a values for the 2SLAQ ensemble alone (offset by a 
factor of log Act = —0.02) and the green stars represent th e 
more strongly clustered 2SLAQ LRGs (|Ross et al.ll2006t ). 
The solid line is the er-projection of the double power-law 
£(r) model which was found to be a good description of the 
2QZ £(r). The relation between w p (a) and £(r) is given by: 



w p (a) 



\f>' 



zdr 



(4) 



The dashed line corresponds to the projection of a 
power law £(r) model, given by £(r) = (r/4.96) -1 ' 85 . 

The fact that the 2SLAQ survey targeted faint QSOs 
is not only an advantage for studies of the luminosity- 
dependence of QSO clustering, but also for 2-space distor- 
tion analyses. The higher spatial density of the combined 
QSO sample should, in principle, improve our statistics 
when studying 2-space distortions, and, in particular, the 
estimation of and P(z) from dynamical and geometri- 
cal £(ct, 7t) distortions. The £(ct, it) measured from the whole 
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Figure 7. The red circles are the w p (cr)/cr measurements for 
the 2QZ+2SLAQ sample. These are very similar to those of the 
2QZ sample alone (blue triangles; da Angela et al., 2005). The 
green stars rep resent the higher clustered LRG sample from the 
2SLAQ survey llRoss et al.l2006l) . The models shown represent the 
projection of a single (dashed line) and a double (solid line) power 
law models. Note that the w p (o)/a values from the individual 
2QZ and 2SLAQ samples have been offset by log Ac of 0.02 and 
—0.02, respectively. 



QSO sample is shown in Fig. [8] (solid contours). The dashed 
lines refer to the 2QZ measurement. 



4 PARAMETER CONSTRAINTS FROM 
REDSHIFT-SPACE DISTORTIONS 

There are basically two mechanisms leading to dynamical 
z-space distortions. As structures grow through gravity, the 
infall of objects to higher-density regions contributes to the 
measured redshifts. If these are assumed to be solely due to 
the Hubble flow, then the large-scale distribution will appear 
flatter, or thinner, along the line of sight, thus "distorting" 
the clustering signal. At smaller scales, the random peculiar 
motions of the objects will also contribute to the measured 
redshifts, and hence distort the measured clustering signal 
for close pairs of objects. If the distribution of distant objects 
has, on average, a spherically symmetric clustering pattern 
in real space, but large velocity dispersion, then the clus- 
tering signal measured in z-space will be smeared along the 
line-of-sight. These features are often referred to as "fingers- 
of-God" , and are commonly seen as elongated structures in 
radial wedge plots of distant galaxy surveys, such as the 
2dFGRS. 

Peculiar velocities are not the only effect leading 
to anisotropies in t he cl ustering pattern. As shown by 
lAlcock fc Paczvnskil ljl979T ). if one assumes a cosmology dif- 



Figure 8. £(cr, it) measured for the 2QZ+2SLAQ sample (solid 
contours) and for the 2QZ sample alone (dashed contours). The 
two measurements are very similar. 



ferent from the true, underlying cosmology of the Universe 
to convert redshifts into distances, the effect on separations 
along the line of sight differs from that affecting the separa- 
tion in the angular coordinate. As a consequence, the clus- 
tering signal might appear elongated (or squashed) in the 
redshift direction. As shown by those authors these geomet- 
ric distortions can be a powerful cosmological test, namely 
to determine fij^. 

Due to their significance at high-z, these potential geo- 
metric distortions have been used to c onstrain cosmologica l 
parameters using QSO catalogues (e.g. lOutram et all l2004); 
21 cm maps of th e epoch of reionisat ion l|Nusserll2005l ) or the 

Lyman a forest (I Becker et al.|[20o3 ). 

However, and as discussed in detail in iBallinger et alj 

(1996), it is sometimes not trivial to disentangle the effects 
of geometric distortions from those caused by peculiar veloc- 
ities. If both the infall parameter /3 and cosmological param- 
eters as or fi° are left as free variables, we expect to see 
a degeneracy between the anisotropies caused by the large 
scale infall and the geometric distortions. Those authors de- 
fine a "flattening factor", which determines, as a function of 
redshift and cosmology, the level of asymmetry expected to 
be seen as a result of geometric distortions, and found that 
its value is degenerate with that of j3. 

The fitting of the dynamical and geometrical distor- 
tions in £(o~, 7r) is de scribed in detail in section 7.7 of 
Ida Angela et al.l (|2005l ). In summary: 

1) for a given value of f3(z), a £((t, 7t) model is gener- 
ated in a chosen test cosmology, through (|Matsubara fc Sutol 
Il996h : 
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fa*) 



(l + |/3W + i/3W 2 )Co(r)Po(/i) (5) 
(|^) + i/3( Z ) 2 )C2(r)P 2 ( M ) (6) 



(7) 



where n is now the cosine of the angle between r and n and 
Pi(l-i) are the Legendre polynomials of order /. £o(V), (,2(r) 
and (,i(r) are the moments of order 0, 2 and 4 of the linear 
£(r) and their form depen ds on the £(r) model adop ted. In 
general, they are given by l|Matsubara fc Sutdfl996h : 



fc»(r) 



(-1)' 

r 2i + l 



da; ) a; 2 ' f — ^ x£ (x) 

I \dx x ) 



(8) 



2) The £(<7, 7r) model is then convolved with the pair- 
wise peculiar velocity distribution to include the small scale 
z-space effects due to the random motions of the QSOs: 



£(o-,tt)= / £{o,-K-w z {l + z)/H{z))f{w z )dw z 



(9) 



where the pairwise veloci ty distribution f(w z ) can be well 
described by a Gaussian |Ratcliff ell 19961 ): 



/K) = 



2tt < iol > 1 /2 



exp 



i Kl s 



2 < wi > 



(10) 



3) Then, the separations a and 7r are scaled to the same 
cosmology that was assumed to measure the actual data. 
The final model for £(cr, 7r) is then compared to the data. 

The relation between the separations a and n in the 
test and assumed cosmologies (referr ed to by the subscript s 
t and a, respectively) is the following (|Ballinger et al.lll996h : 

f* = f± cr a = "f^CTa (11) 

A* 

7Ti = /||1To = -T-7T (12) 

where ^4 and B are defined as (for spatially flat cosmologies) : 

(13) 



.4 



dz 



H J ^n A + n° m (i + z>y 



(14) 



In the linear regime, the correlation function in the as- 
sumed cosmology will be the same as the correlation function 
in the test cosmology, given that the separations are scaled 
appropriately, i.e.: 



£t(Ct) 1"t) = £a(o"a, TTa)- 



(15) 



4) This method is repeated for different test cosmolo- 
gies and values of /3(z). 

Given the similarities between the £ (a, n) contours in 
Fig. [8] in addition to very similar £(s) and w p (a) mea- 
surements, we would not expect the constraints put on 
(3(z) and from the 2QZ+2SLAQ dynamical distortions 
to differ from those ob tained from the 2QZ sample alone 
(|da Angela et al.ll2005h . assuming that all the underlying 
assumptions remain the same (e.g., £(r) shape and ampli- 
tude, velocity dispersion, scale-independent bias). We now 



repeat the method adopted for fitting the 2QZ dynami- 
cal and geometrical distortions, but also utilising the new 
2SLAQ ensemble. The question now arises if the same 
£(r) model should be assumed, or if the velocity disper- 
sion of the QSOs should still be fixed at 800 km s -1 . It 
can be seen in Figs. [S] and [7] that the 2QZ double power- 
law £(r) model is still a good description of both £(s) 
and w p (<j) measurements for the combined sample. As the 
2QZ and 2SLAQ samples have similar N(z) we would not 
expect to see clustering differences between them due to 
redshift evolution. Any potential clustering difference be- 
tween both sets would be due to the different luminosity 
of the sa mples. However, as suggested by both observa- 
tions (e.g. ICroom et all 120051; lAdelberger fc Steidell [20051 : 



IPorciani fc Norberd 12006 : I Myers et aTT 20061 ). andsimula- 
tions |Lidz et alj|2006t iHopkins et al.ll2006l ) and, more im- 
portantly, as we shall see later in this paper, QSO clustering 
is very weakly luminosity-dependent. We therefore assume 
the same double power-law £(r) prescription as used for the 
2QZ sample. We also assume the same velocity dispersion as 
for the 2QZ sample alone. It is not unlikely that the 2SLAQ 
QSOs would have, o n average, a dif f erent ve locity dispersion. 
As po inte d out by Berlind et all (|2003l ) . lYoshikawa et al] 
l|2003l) , or iTinker et all (|2006l ). galaxies can be a biased 
tracer of the dark matter velocity distribution, just as they 
are of the dark matter spatial distribution. However, as 
found for the 2dFGRS galaxies a nd predicted by HO D (Halo 
occupation distribution) models l| Tinker et alJliooil ) , the ex- 
pected difference for Mtj^— 20 is not significant. In addi- 
tion, as most of the 2-error is due to measurement erro r 
rather than intrinsic velocity dispersion (|Croom et al.ll2005l ). 
we chose to continue assuming < w\ y 1 / 2 — 800kms _1 . 

The fit to the distortions in £ (a, n) was performed with 
the same assumptions and over the same range of scales as in 
the previous 2QZ analysis. The result is shown in Fig. [9] As 
expected, the contours are indeed tighter than the ones ob- 
tained when fitting only the 2QZ £(<r, ir). This is due to the 
increased number of pairs, not only from the 2SLAQ sam- 
ple alone but also from the cross-pairs in the NGC between 
the two ensembles, as they probe overlapping volumes. Also 
shown are the la and 2a confidence levels predicted from 
clustering evolution and linear theory of density perturba- 
tions (dashed lines). The dotted line is, as usual, the la joint 
confidence levels from both constraints. The best fitting val- 
ues are Q.^ = 0.25±g:{£, (3(z) = 0.60^;", corresponding 
to a Xmin = 1-02 (12 d.o.f.). Although these results favour 
a somewhat higher value of /3 than the previous 2QZ only 
£(<r, 7r) constraint, both obtained results are self-consistent, 
within the associated errors. We should point out that the 
size of the error bars does not take into account any po- 
tential correlation between £(a, n) bins but this is expected 
to be small. Finally the above derived values for fi^ and 
(3 imply a value of the QSO bias of b(z = 1.4) = 1.5 ± 0.2 
which is slightly lower than, but not inconsistent with, the 
value of biz — 1.4) « 2 derived below from purely the QSO 
clustering amplitude. 



5 LUMINOSITY-REDSHIFT DEGENERACY 

A few recent wo rks have looked at t he evolution o f QSO 
clustering (e.g. ICroom et ail 120051 : IPorciani etafl 120041 : 
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Figure 9. Confidence levels in the [O^, (3(z)] plane obtained from 
fitting the 2QZ and 2SLAQ £(cr, n) z-space distortions (solid lines 
and shaded contours). Dashed lines show the 1 a and 2 a con- 
straints from linear theory evolution. The dotted contour is the 
1 a joint confidence level. 



I Myers et alj l2006t ). These suggest an increase of QSO 
clustering amplitude with redshift, a trend which is more 
significant at z> 1.6. This evolution contrasts with that 
expected from a long-lived QSO population model, or linear 
theory predictions, which generally predict a decrease of 
clust e ring a mplitude with increasing redshift l|Croom et al.l 
l200ll . 120051 ). The range of magnitudes covered by the QSO 
surveys used in these studies has not fully permitted the 
study of the luminosity dependence of QSO clustering. 
However, the combination of the 2QZ and 2SLAQ samples 
probably sees its greatest scientific contribution precisely 
in the range of luminosities it probes and for the first time 
allows a more rigorous dete rmination of t h e QSO clustering 
dependence on luminosity. ICroom et al.l (2002) have used 
the 2QZ sample alone for this purpose. Their results suggest 
that additional, fainter data, such as those obtained for 
2SLAQ, are essential to pursue this goal. 

To estimate the bj band absolute magnitude, Mb,, we 
compute: 



Mb j (z) =bj- Kb j (z) - A bj + 5-5 log(d), 



(16) 



where bj is the apparent magnitude, K bj the k-correction in 
the bj magnitude, Ab, the dust correction and d the lumi- 
nosity distance that corresponds to the redshift z, measured 
in parsecs. The value of the k-correction was taken from 
ICristiani fe Viol (ll990T l. The galactic dust cor rection, Ab , is 
deter mined through: A bj = A.035E(B - V) (|Schlegel et all 
1 19981 ). 

The above formula is used to determine the absolute 
magnitude of the 2QZ QSOs. To include the dust correction 
when determining the absolute magnitude of the 2SLAQ 



QSOs, one subtracts the g magnitude galactic extinction 
(ffred) at the QSO's coordinates from the observed apparent 
magnitude (g): g' — g — g rec l> where g' is the dust-corrected 
<?-band QSO magnitude. The other subtlety in combining 
the two QSO samples is accounting for the relation between 
the observed bj and g magnitudes. However, this becomes 
quite simple as the transmissivity curves of the filters have 
a significant overlap and the same zer o-point. Thus, we ca n 
treat these bands as being equivalent (|Richards et aljfeoosl ) . 
Hereafter, and for the sake of simplicity, we shall refer to the 
QSO absolute magnitudes for both samples as if they had 
been measured in the bj band, and represent both of them 
as Mbj. Therefore, the 2SLAQ QSOs' absolute magnitude 
is determined by: 



M bj (z)&g'-K bj (z)+5-5 log(rf) , 



(17) 



where g' already includes the dust correction in the g band. 

Fig. [TU] shows how the 2QZ and 2SLAQ are distributed 
in the [Mbj, z] plane. The 2QZ QSOs are shown in red and 
the 2SLAQ in blue. The cyan lines represent the adopted 
2QZ bj < 20.85 and 2SLAQ 20.5 < g < 21.85 magnitude 
cuts. The QSO samples span the z-range 0.3 < z < 2.9. 
The yellow line shows how M bj changes with 2. We adopted 
a second-order polynomial model to determine M£ (z ) 
l|Bovle et al.ll2000l : ICroom et ai]|2004l : iRichards et~ail2005l ) : 

Mbj (z) = Mbj (0) - 2.5(fci2 - k 2 z 2 ) (18) 

We adopt the values obtained by ICroom et all (12004): 



Mb , (0 ) = -21.61, fci = 1.39, k 2 = -0.29. IRichards et al.l 



(2005) showed that the parameterisation of the M^(z) 
model is only marginally affected by including or not the 
2SLAQ QSOs. The yellow line in Fig. \W\ only extends to 
z = 2.2 given the fitting range used in this parameterisa- 
tion. 

The flux-limited nature of these two surveys is evident 
in this plot. More luminous QSOs lie at higher redshifts 
while fainter ones have lower redshifts. This means that, 
unless we probe a wide window in magnitude-space with 
our QSO surveys, it will be intrinsically hard to determine 
how QSO physical properties change with luminosity, for a 
fixed redshift. By combining the 2SLAQ and 2QZ samples 
we are widening the magnitude window and hence making 
it possible to determine the dependence of QSO clustering 
on luminosity, free of any evolutionary effects. 

Fig. [TD] also shows how, using the two surveys together, 
we can look at a specific redshift range and determine the 
QSO clustering in different magnitude samples. This "verti- 
cal approach" to the [Mbj,z] distribution is possibly more 
physically justifiable than simply analysing QSO cluster- 
ing dependence on redshift or apparent magnitude. Tests of 
models where comparisons at fixed luminosity are required 
certainly need as full coverage as possible of the luminosity- 
redshift plane. Indeed, the long-lived QSO model has been 
easiest to test in previous samples, since comparing intrin- 
sically low luminosity QSOs at low redshift with high lu- 
minosity QSOs at high redshift makes more sense in a PLE 
model where the two are hypothesised to be directly related. 
These results have been used to argue against a long-lived 
model for QSOs with the 2QZ results of £20 (see equn. [T9l 
below) appearing to rise, if anything, r ather than fall with 
increasing redshift (|Croom et al.ll2004 ). However, the low 
redshift (z « 0.02) IRAS selected Seyfert 1 and 2 results of 
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Figure 10. Magnitude and redshift bins adopted for 2QZ and 
2SLAQ QSOs. The numbers in each division of the "grid" are 
the numbers of 2QZ and 2SLAQ QSOs in the specific z and M b 7 



Figure 11. QSO £(s) measured in different magnitude and red- 
shift bins. The order of the panels is the same as that of the 
[Mj, , z] intervals in Fig. 1101 The dashed line shows the best fit- 
ting power-law to the £(s) of the full sample. The solid line is the 
£(s) power-law fit to the data in each individual panel. 



iGeorgantopoulos fc Shanksl i|l994 ) give £>o = 0.52 ± 0.13 
in goo d agreement with the low redshift SDSS AGN re- 
sults of I Wake et all (|2004r ) which give £20 = 0.48 ± 0.03 and 
adding these r elativ ely high amplitude points to Fig. 21a of 
ICroom et al.l (|2004 ) may make the clustering case against 
the long-lived model less strong. We note also that the low 
value of £20 = 0.22 ±0.08 measured for the cross-correlation 
of low luminosity QSOs with Ly man-break galaxies at z=2.5 
by I Adelberger fc Ste idcl (2005]) also goes against the trend 
for higher clustering amplitudes for higher redshift QSOs. 
But whatever the hypothesis for QSO lifetime, the extended 
luminosity range of the 2SLAQ sample means that we can 
now test the generic prediction of these 'high peaks' bias 
models for higher clustering amplitudes for more luminous, 
rare QSOs at fixed redshift. 



6 CLUSTERING AS A FUNCTION OF 
MAGNITUDE AND REDSHIFT 

Dividing up the QSO samples into magnitude and red- 
shift bins significantly increases the error on our cluster- 
ing measurements, simply due to the much smaller num- 
ber of objects in each bin compared to the total number 
of QSOs (numbers in Fig. I10p . This is also evident in Fig. 
111! where we plot the £(s) measurements in each of the 
panels in Fig. 1101 The dashed line shows the best fitting 
power-law model to the overall 2QZ+2SLAQ sample, over 
the 3 < s < 50 7i _1 Mpc range (£(s) = (s/6.20)" 1 ' 66 ). The 
solid lines are the best power-law models to each individual 
[Mbj,z] interval, fixing the £(s) slope to 7 = 1.66 and per- 



forming a x 2 fit to determine the amplitude. The order of 
the panels in Fig. [11] is the same as in the panels presented 
in the [Mb , , z] plane in Fig. \W\ 

By visually comparing the dashed and solid lines, we 
observe no dependence of QSO clustering on luminosity nor 
redshift. However, the size of the errorbars motivates the fur- 
ther use of more statistically robust tools. We therefore use 
the integrated correlation function up to 20 h~ Mpc in order 
to quantify the clustering amplitude in each magnitude-z 
bin. This quantity is then normalised to the volume con- 
tained in a 20 /i -1 Mpc sphere: 

60 = J tt s > 2ds ( 19 ) 

The choice of using 20 ft" 1 Mpc as the radius of the 
spheres to compute the averaged correlation function is due 
to the fact that this is a large enoug h scale for lin e ar the - 
ory to be applied and, as shown by ICroom et al.l ([2005), 
small-scale z-space distortions do not significantly affect the 
clustering measurements, when averaged over this range of 
scales. In addition, and as seen in Fig. |6j we can estimate 
the uncertainty through computing Poisson errors, and scale 
this by a factor of 1.25. This estimate should provide a fair 
description of the uncertainty on the correlation function 
measurements, and significantly reduce the computing time. 

We computed £20 using the Hamilton estimator in each 
of the bins shown in Fig. 1101 The results for each redshift 
slice are shown in the four panels in Fig. 1121 Red circles show 
the measurements in each magnitude bin. The shaded grey 
area shows the la £20 measurement for QSOs of all luminosi- 



The 2SLAQ survey: QSO clustering and the L-z degeneracy 11 



3 2 10 -1-2-3-4-5 





1,2 
1.0 
0.8 
0.6 
0.4 
0.2 



4 



-20 -22 -24 -26 -28 -30 



-20 -22 -24 -26 -28 -30 



Figure 12. The four panels represent the §20 measurements in 
different redshift bins. The median redshift of each z-interval is 
indicated in the top left of each graph. The top horizontal axis 
shows the absolute magnitude difference, relative to M£ (< z >). 
The red circles are the §20 measurements in different absolute 
magnitude bins, and are centred on the median values of each 
bin. The shaded area is the lcr £20 interval for all the QSOs in 
that specific redshift interval. The horizontal length of the shaded 
area represents the range of Mbj values of QSOs in that redshift 
interval. 



ties in that specific redshift slice and its length indicates the 
total range of magnitudes included. The dashed line repre- 
sents the average value of £20, for all redshift and magnitude 
ranges. It should be pointed out that the bin sizes were cho- 
sen in such a way that the precision of the clustering mea- 
surements was maximised, and therefore the distribution of 
QSOs in a given z-slice is not constant for all magnitudes. 
Thus, we do not expect our £20 measurements to be equidis- 
tant along the horizontal axis, as these are centred on the 
median values in M b] of each bin. The top axis indicates 
the magnitude difference with respect to M bj (< z >), at 
the median redshift of that specific "z - slice" . The "rising" 
of the grey area as we move to h i gher r edshifts is consistent 
with the results of ICroom et a l. (2005), who also found an 
increase of clustering amplitude with redshift, for the 2QZ 
QSOs. 

The number of QSOs in each M bj — z bin, indicated in 
Fig. 1101 is now reflected in the sizes of the £20 error bars. In 
the first, lower-z panel, for instance, the Mbj — z bin with 
only 533 QSOs corresponds to the £20 measurement with 
the largest error bar. The two intermediate z-slices are the 
ones where most of the gain of the 2SLAQ is observed and 
the ones with highest statistical value. Our results are in 
agreement with the hypothesis of a luminosity-independent 



QSO clustering being constant with redshift and luminosity 
is not supported by the data (Xred = 2.50). 



7 BIAS AND HALO MASSES 

The £20 vs. Mb j results motivate the analys is of the de- 
pende nce of bias on luminosity and redshift. ICroom et al.l 
(2005) investigated the redshift evolution of QSO bias, us- 
ing the 2QZ survey data. They found that the QSO bias 
does evolve very strongly with redshift; as the mass cluster- 
ing amplitude decreases with increasing redshift, the slight 
upward trend observed in the 2QZ £20 reveals a strong in- 
crease of bias with z. 

Under the assumption of a sc ale- independen t bias, the 
bias can be obtained through (e.g. |Peebleslll98(ih : 



Mr) 



fofr,20) 
0,20)' 



(20) 



where £q(r, 20) and £p( r > 20) represent the QSO and mat- 
ter real-space correlation functions, respectively, averaged in 
20/i _1 Mpc spheres. The z -space and re al-space correlation 
functions can be given by (lKaiserj ll987Tl: 



Cq(s,20) 



-/?+ V 
3 M 5 M 



£«(»■, 20) 



(21) 



Combining both equations and taking into account that 
P — Q, 1 ^ 6 /b leaves us with a quadratic equation in b. Solving 
it (see ICroom et al.ll2005T ) leads to: 



b(z) = 



£q(s,20) 4^ 2 (z) Q° m b (z) 
£p(r,20) 45 3 



(22) 



clustering {\ r 



1.16, over 12 d.o.f). The hypothesis of 



Therefore, we can use our £q(s, 20) measurements, rep- 
resented in Fig. I12l and. together with a theoretical estimate 
of £p( r > 20), determine the bias that corresponds to that the- 
oretical assumption and the observed clustering measure- 
ments, on the assumption of a cosmological model. Our re- 
sults are shown in Fig. 1131 T o estimate £ p ( r, 20) , we use the 
P(k) non-linear estimate of ISmith et al.l ([2003). To deter- 
mine £p( r ) we Fourier transform this P(k) estimate, and in- 
tegrate the result up to s ^ 20 /i _1 Mpc to compute £ p (r, 20). 
The parameters used to generate the P(k) model were: 
fim = 0.3, fi° = 0.7, r = 0.17 and, for a better compar- 
ison with Croom et al.'s (2005) results, as = 0.84. This 
value is consistent with recent studies (e.g. iPercival et al.l 
l2002l : lTvtle~ et al. I [2004m . even though recent measurements 
also tend to suggest somewhat lower values (|Spergel et al.l 
l200fj) . 

The stars in Fig. [13] represent the b estimates for the 
magnitude-integrated samples, corresponding to the shaded 
areas in Fig. 1121 The se values are v ery m uch in agreement 
with those found by ICroom et alj (|2005l ). using a similar 
method. The dashed line is the empirical description of 

b(z) = 0.53 + 0.289(1 + zf 

found by those authors. 

The circles refer to our measurements in different mag- 
nitude bins. The red ones correspond to the faintest, Mbj > 
-23.5 QSOs; the blue ones to the -24.5 < M bj < -23.5 
range; the green circles represent the QSOs with —25.5 < 
M bj < -24.5 and the brightest, M bj < -25.5 QSOs are 
represented by the yellow circles. Given the size of the error 
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Figure 13. Bias evolution for different luminosity QSOs. The 
different colours refer to different absolute magnitude bins. The 
stars are the result for all the QSOs in each specific rcdshift bin. 
The dashed line is the empirical model of Croom et al. (2005). 
Each point is represented in the median redshift of all the QSOs 
in the specific Mj,, and z-ranges. 



bars, which are related to the errors on the associated £(20) 
measurements, no categorical conclusion can be drawn, re- 
garding the possibility of a luminosity-dependent QSO bias. 
The uprise in the bias values with redshift is unrelated to 
the different QSO luminosities, as a somewhat positive trend 
occurs for all QSOs irrespective of their magnitude. This is 
not entirely true for the brightest, Mbj < —25.5, QSOs, 
(in yellow) for which the bias at z ~ 1.3 seems higher 
than at higher redshifts. However, given the small number 
of QSOs (355) within that redshift/magnitude range, this 
result would need further study. 

The b values for each magnitude are centred in the me- 
dian redshift of the QSO sub-sample from which b was deter- 
mined. Hence, in each redshift bin, the 2-displacement of dif- 
ferent magnitude points is due to the non-uniform distribu- 
tion of the QSOs in the [Mbj, z] plane. That z-displacement, 
together with the colour-code on the left side of the plot, 
makes it easy to relate Fig. [12] to Fig. 1131 

The red, 2SLAQ-dominated, fainter bin at z ~ 1 with 
a relatively small err or bar deviates sig nificantly from 
the empirical model of ICroom et al.l (|2005l ). However, the 
overall trend is conservatively consistent with a luminosity- 
independent QSO bias. 

The bias of the QSOs is related to the mass of the 
dark matter halo they inhabit. In a Gaussian random field 
the higher the fluctuation threshold the higher the clus- 
tering amplitude of fluctuations. Therefore, by measuring 
the clustering of QSOs we can infer the mass of the haloes 
the QSOs inhabit. The formalism relating bias and halo 



mass was firstly developed bv lMo fc White) (|l996l ). who as- 
sumed a spherical collapse model. This was then extended 
to more compli c ated geometries, such as ellipsoidal collapse, 
bv lSheth et al.l (1200 ll ). In the analysis in this work the latter 
will be the adopted formalism. According to these authors, 
the bias can be related to the dark halo mass by: 



6(Mi 



DMH, Z) 



1 + 



\/~a(ai> 2 ) + \fab{av 2 ) 1 



(a!/ 2 ) c + fc(l-c)(l-c/2) 



(23) 



with a = 0.707, 6 = 0.5 and c = 0.6. v is denned as v = 
8 c (z)/a(MnMH, z). S c is the critical density fo r collapse, and 
is giv en by: S c = 0.15(127r) 2/3 ft m (,z) - 0055 (|Navarro et all 
1 1997ft . <t(M D mh, z) = a(M DMH )G{z), where a(M UMH ) is the 
rms fluctuation of the density field on the mass sc ale with 
value tfpMH and G(z) i s the linear growth factor dPeeblesI 
ll984ICarroll et al.| [l992). <j(M D mh) can hence be computed 



<x(M D 



1 

2^ 



k 2 P(k)w(kr) 2 dk 



(24) 



where P(k) is the power spectrum of density perturbations 
and w(kr) is the Fouri er transform of a spherical top hat, 
which can be given by |Peebleslll980ft : 



w(kr) = 3 



sin(fcr) — krcos(kr) 



(kr) 3 

where the radius r is related to the mass by: 
r= f 3A/dmi X ' " 

I 47T/90 



(25) 



(26) 



and po = SlmPci it i s the present mean density of the Uni- 
verse, given by p = 2.78 x lO 11 fi£ 1 h 2 M Mpc -3 . 

Here, we adopt a linear form of the power spectrum, 
P(k) = PoT(k) 2 k n , where Po is simply a normalisation pa- 
rameter that depends on <rg and T(k) is the transfer func- 
tion, which we descr ibe through the analytical formula of 
iBardeen et all l| 1986ft . 

The results of performing this analysis using our deter- 
mination of the bias is shown in Fig. Q3] The panels show 
the dark matter halo mass associated with different luminos- 
ity QSOs, in the same redshift intervals as those plotted in 
Fig. 1121 The horizontal axes show the QSO absolute magni- 
tude (bottom), and its difference relative to Afj (top axis), 
similarly to Fig. 1121 In each panel, the red circles represent 
the Mdmh measurements in different magnitude bins, with 
error bars being the uncertainties corresponding to those 
obtained in our previous b(z) estimates. The increase in the 
relative errors in Fig. Q3] as compared to Figs. 1121 1131 is due 
to the relatively flat slope of the cf(Mdmh) relation for the 
AC DM model as obtained from equn. [24l The shaded areas 
represent the la Mdmh confidence levels when estimating 
the masses associated with all QSOs, irrespective of their 
luminosities. 

We find that, at all redshifts, QSOs seem to inhabit 
Mdmh ~ 3 x 10 12 Mq haloes (das hed line), very muc h 
in agreement with what was found bv lCroom etai] l|2005ft . 
As pointed out by those authors, this result appears to dis- 
favour the picture of a long-lived QSO population. As the 
dark matter halo masses grow, with decreasing redshift, we 
would expect to see lower-z QSOs in more massive haloes, if 
that were the case. The fact that we do not, means that at 
consecutive redshift intervals, we may not be observing the 
same QSO population, but rather distinct sets of objects. 
However, this conclusion is based on the 2QZ and 2SLAQ 
results alone and so the caveat made at the end of Section 
5 about the higher clustering amplitudes measured by other 
authors for low redshift AGN still applies. 

We also find, through our results, no evidence for Mdmh 
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Figure 14. The four panels show the Mdmh estimates in differ- 
ent redshift bins. The median redshift of each ^-interval is indi- 
cated in the top left of each graph. The top horizontal axis shows 
the magnitude difference relative to Mi (< z >). The red circles 
show the dark matter halo mass measurements in different abso- 
lute magnitude bins, and are centred on the median values of each 
bin. The shaded area is the la interval for the Mdmh value of all 
QSOs in that specific redshift interval. The horizontal length of 
the shaded area represents the range of Mj, 7 values for the QSOs 
in the redshift interval. The dashed line shows the average AfrjMH 
at all redshifts. 



segregation with QSO magnitude at fixed redshift. All the 
values seem to be consistent with a flat Mdmh — Mbj trend, 
indicating that QSOs seem to live in ~ 1O 12 /i _1 M haloes, 
independently of their luminosity. This behaviour is incon- 
sistent with simple, 'high peaks', models of QSO biasing 
where rare, luminous QSOs might be expected to occupy 
higher mass haloes. 



with different luminosity QSOs, given that we know the mass 
of the haloes that they inhabit, and thus determine if indeed 
more luminous QSOs are associated with more massive black 
holes. For each of these two e volutionary sc e narios , and fol- 
lowing iFerraresei (|2002h and ICroom et all (|2005h . we will 
consider three possibilities for the dark matter halo profile, 
which affect each of assumed scenarios differently. We will 
consider: a) an isoth ermal dark matter profile; b) a NFW 
(|Navarro et al.l fl997T ) profile and c ) a profile inferred from 
weak lensing studies l|Seliakl [2002). which, for the sake of 
simplicity, we will refer to as the "lensing" profile. 

When assuming a ^-independent A/bh — Mdmh cor- 
relation, the three possible (a), b) and c)) halo p rofiles 
correspond to the following relations (jFerraresd 12002 ) : 



1. a) Isothermal profile: 



Mb 



0.027 



M D 



lO 8 Af V io 12 m q 

1. b) NFW profile: 

Mbh I Mdmh 

1O 8 M ~ ^10 12 M© 

1. c) "Lensing" profile 
Mbh „ „„ / A'/dmh 



1O 8 M 



0.67 



10 12 A/ ( 



(27) 



(28) 



(29) 



If we assume a z-independent correlation between the 
black h ole mass and the c ircular velocity in the associated 
bulges (IShields et alj|2003t 2.) then other relations are ob- 
tained. Following Croom et alj (|2005l ) and IWvithe fc Loebl 
(2005aj), the equivalent relations between the dark matter 
halo mass and the black hole mass are given by: 



Mbh = e 



M D 



2/3 



1O 12 M 
where A has the form: 
A = 18tt 2 + 82 (n ra (z). 



AD." 



5/6 



187r 2 a, 



39 (Q m (z) 



(1 



5/2 



(30) 



(31) 



The constant e is related to the halo density profile. 
Different values of e will correspond to the s ame scenarios as 
conside red in case 1.. Hence, and following IWvithe fc Loebl 
|2005al ). we have that: 



8 ESTIMATING BLACK HOLE MASSES FOR 
DIFFERENT LUMINOSITY QSOS 

Several models and theoretical studies have been developed 
to try to determine the relation between the mass of the 
dark matter halo and the mass of the black holes associated 
with the observed QSOs. Here we will consi der the two pos- 
sible ev olutionary scenarios considere d by Wyithc fc Loebl 
(2005a), both based on the results of IFerraresei (2002) 1 : 1. 
a correlation exists between the dark m atter halo mass 
(Mdmh) and the black hole mass (M BH ) (|Ferraresel [20021 
and this relation is unevolving with redshift; 2. instead, the 
correlation between the bulge velo city dispersion (or circular 
velocity) and the blac k hole mass (|Ferrarese fc MerrittfeOOOl ; 
iGebhardt et~ai]|2000l ) is assumed to be unevolving with red- 
shift. We can then estimate the black hole masses associated 



2. a) For an isothermal profile: 

1Q -5.1 

2. b) For a NFW profile: 
3.7 x 10" 51 

2. c) For the "lensing" profile: 
25 x 10~ 5 1 



(32) 



(33) 



(34) 



Again, as in case 1., the three different possibilities con- 
sidered for the density profile differ only in terms of a nor- 
malisation parameter, in this case, given by the constant 
e. 

We now use relations 1. - 2., a), b) and c), to determine 
the mass of the black holes that correspond to our Mdmh 
measurements, under different assumptions, and determine 
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Figure 15. Black hole mass as a function of luminosity, in differ- 
ent redshift bins. The filled symbols and solid lines are obtained 
assuming a A/bh - A/dmh relation which is independent of z. 
The dashed lines and open symbols, which also correspond to the 
errorbars with larger tickmarks, assume a z-independent A/bh - 
cT c relation. In both cases, the circles, squares and triangles cor- 
respond to isothermal, NFW and lens-studies-based halo density 
profile, respectively. The points are located at the median lumi- 
nosity value of the QSO sub-sample to which they correspond. 
On the vertical axis on the right of each panel is the equivalent 
Eddington luminosity scale to that for A/bHi on the left. The yel- 
low shaded area represents the super-Eddington, L/LEdd > lj 
regime. The dashed yellow line corresponds to a Eddington effi- 
ciency e = 0.01. It can be seen that some models imply super- 
Eddington solutions, and hence are unlikely to occur. Most of the 
models though, correspond to 0.01^ e<> 1.0 values. 



if, with the current data, we can relate the black hole mass 
to the QSO luminosity. 

Our results are shown in Fig. 1151 Each panel shows the 
results obtained in a given redshift bin. Plotted is the black 
hole mass as a function of QSO luminosity. To determine 
the b olometric luminosity from we use l|Croom et al.l 
|2005| V. 



_ 10 (79.42-Af bj )/2 



•w 



(35) 



The blue filled symbols and solid lines refer to hypoth- 
esis 1., where we assume a A/bh — Mdmh z-independent 
relation. The red open symbols and dashed lines relate to 
hypothesis 2., where we assume a Mbh — cr c relation inde- 
pendent of z. The filled and open circles show the a) esti- 
mates, in Eqs. 1271 and 1321 respectively, on which we assume 
an isothermal density profile. The squares show the results if 
we assume a NFW profile (b)) and the triangles if we assume 
the lensing profile (c)). The error bars are the corresponding 
uncertainties to those on the Mdmh measurements, plotted 
in Fig. 1141 To distinguish between the error bars, the ones 



that refer to hypothesis 1. are represented with short tick 
marks, whereas the ones that refer to hypothesis 2. have 
longer tick marks. 

For both of the assumptions, 1. or 2., the dark matter 
halo "lensing" density profile corresponds to more massive 
black holes, and the isothermal density profile corresponds 
to the least massive black holes, as expected. Also, it be- 
comes evident that assuming different profiles, being it un- 
der z-independent A/bh — Mdmh or A/bh — o"c scenarios, 
simply "shifts" the A/bh — log(L) relation vertically. Even 
though the errors associated with the Mbh are large, we can 
say th at our values are consistent with those of lCroom et al.l 
l|2005h . who studied the evolution of Mbh with redshift. 

Also shown, on the right-hand side vertical axis in each 
panel, is the Eddington luminosity. This is determined di- 
rectly from the black hole mass as follows: 



Ll 



irr 



Mb 



10 8 M Q 



W 



(36) 



The yellow area in the bottom of each panel represents 
the values of A/bh that correspond to "super-Eddington" 
solutions ie, L/LEdd > T The dashed line represents 
the Mbh — log(i) relation for an Eddington efficiency of 
e = L/LEdd = 0.01. It can be seen that, for some of the sce- 
narios considered, the mean efficiency is super-Eddington, 
in particular for models l.a) and l.b), ie, assuming an 
isothermal profile and an NFW profile, when considering 
that the Mbh — Mdmh relation that does not evolve with 
redshift. One could argue that these relations are therefore 
unlikely to occur. Most of the remaining models suggest 
accretion efficiencies of 0.01< e< 1. It is somewhat unfor- 
tunate that the size of error bars do not allow us to draw 
conclusions regarding the significance of potential changes 
of black hole mass with luminosity of the associated QSO. 

We averaged the data over the whole redshift range to 
test, through a simple \ 2 analysis, the hypothesis that QSOs 
do not accrete at a fixed fraction of Eddington. Fig. [16] rep- 
resents the results, by assuming the "lensing" halo density 
profile and z-independent Mbh - cr c relation (open red tri- 
angles). Also shown is the best fitting value of e for that 
assumption. 

From the "flat" trend observed in the measured A/bh— L 
relation, black hole mass seems approximately independent 
of QSO luminosity. However, this does not permit us to ex- 
clude the hypothesis that high-z QSOs accrete at a fixed 
fraction of Eddington, as a model characterised by a con- 
stant value of e is still a good fit to the data ( e = 0.044 with 
\ 2 = 1.58; 3 d.f.). Gi v en re cent studies of iHopkins et al.l 
l|2005h and iLidz etafl l^OOrJ ), who argue that bright and 
faint QSOs are similar sources, but observed at different 
stages of their activity, one could expect both luminous 
and faint QSOs to be associated with equally massive black 
holes. This would thus lead to higher values of accretion ef- 
ficiency for brighter QSOs and lower for fainter QSOs. Such 
a model can still be in agreement with the current analysis, 
given the "flat" trend of the Mbh values as a function of 
luminosity. 

Hence, our results show that, if halo mass and black 
hole mass are closely correlated, then we cannot reject a 
model where black hole mass depends on QSO luminosity 
and accretion efficiency. It should be noted that we have as- 
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Figure 16. Black hole mass as a function of luminosity, over all 
redshifts. Here we assume the "lensing" halo density profile and 
that the Mbh - sigma relation is z-independcnt. The best fitting 
value of e is shown by the dashed yellow line. 



sumed that the dispersion in the black hole mass and halo 
mass is small. This assumption is supported by the existence 
of rea sonably tight bulge mass - velocity dispersion rela- 
tions (|Tremaine et al]|2002l ). But clearly, if this assumption 
proved incorrect, then the results above would be affected 
by the high dispersion in the Mb u i ge — Mbh relation. 



9 CONCLUSIONS 

The 2SLAQ QSO survey is an important tool for QSO clus- 
tering studies. Firstly, the 2SLAQ QSO survey complements 
the previous 2QZ sample in terms of 2-space distortion anal- 
yses. We have shown that a double-power law £(r) model, 
which is a good description of the 2QZ real-space cluster- 
ing, still describes well both the z-space and projected clus- 
tering measurements of the 2QZ and 2SLAQ samples com- 
bined. We fit the dynamical and geometrical distortions of 
the £(cr, 7r ) cont ours, extending the formalism developed by 
iHamiltonl (|l992T ) and lMatsubara fc Sutcl (| 19961 ) to include a 
double-po wer law £(r) model and fitt i ng different "test" cos- 
mologies dAlcock fc Paczvnskl Il979l ; iBallinger et all Il99rj ; 
Ida Angela et al.ll2005f ). We find that the subsequent confi- 
dence levels obtained in flj^ and (3(z) are similar to those ob- 
tained when using solely the 2QZ data, but tighter due to the 
increased statistics from extra 2SLAQ QSO pairs, and also 
the additional cross-correlation pairs in the NGC 2SLAQ 
and 2QZ overlapping volumes. When combining these re- 
sults with orthogonal contours obtained from linear the- 
ory of density perturbations, we find that = 0.251q g7> 
(3(z) = 0.60^q '\ 4 L , similar to the values obtained from the 



2QZ data alone (Q° m = 0.35 
new results imply b(z = 1.4) = 1.5 ± 0.2 for the QSO bias. 

Secondly, the 2SLAQ QSO survey constitutes a new 
dataset with a potentially central role in terms of breaking 
the L-z degeneracy. The sample extends 1 magnitude fainter 
than the 2QZ, and spans the same z-range. Hence, the com- 
bination of both provides a unique dataset, as the overall 
magnitude range probed is similar, both at low and high- 
z. This allows us to interpret clustering results and possi- 
ble luminosity-dependent measurements in different redshift 
bins, hence reducing any evolutionary biases. Our results are 
consistent with luminosity -independent QSO c lustering and 
in agreement with those of ICroom et all f 20051 ); QSOs seem 
to inhabit ~3x 10 12 /i -1 Mq haloes, independently of their 
redshift or luminosity. Our results do not show a tight cor- 
relation between halo mass and QSO luminosity at fixed 
redshift, as would be expected from simple "high peaks" 
models of QSO biasing where fainter QSOs populate lower 
mass haloes. 

Our Mdmh vs. Mb j results agree with the predictions 
of lLidz et al.l J2006I). wh o se sim ula tion r esults based on the 
models of iHopkins et all (|2005l laHbl 120061 ) suggest that QSO 
luminosity may not be correlated with the mass of the host 
dark matter halo. The reason is that the same massive haloes 
host faint and bright QSOs and the difference in luminosity 
is due to the QSOs being observed in different periods of 
their lifetime. Another consequence is that QSO clustering 
should not correlate strongly with luminosity, again, just as 
shown by our data. 

These authors' analysis also support the res ults shown 
in the present paper and bv lCroom et al.l ((2005) , and con- 
clude that QSO clustering and halo mass do not evolve 
strongly with redshift, even though QSO bias substantially 
increases as we move to higher z. This could hint at possible 
anti-hierarchical QSO formation (|Merlon i 2005; Cowie et al.l 
120031 ; iLidz et al.ll2006l ). as haloes harbouring QSOs would 
have deeper potential wells at high-z than at low-z, leading 
to more luminous black holes being observed at high-z than 
at low-z. The reason for the rapid decrease of QSO bias with 
time is related to haloes of ~ 10 12 — 10 13 Mq corresponding 
to rarer, high-density-contrast peaks at higher redshift. The 
results of those authors also predict that a large range in 
QSO luminosity should correspond to a very restricted range 
in QSO halo masses, as our observations and measurements 
seem to indicate. 

By assuming different density profiles for the dark mat- 
ter halo and z-independent relations (such as Mbh - Mdmh 
or Mbh - er c ) we can estimate the masses of the black holes 
associated with the QSOs. If the Eddington limit is a rel- 
evant limit for the accretion rate, and if one assumes that 
the Mbh - Mdmh relation is z-independent, then isothermal 
and NFW density profiles are not likely to be appropriate 
for the haloes these QSOs inhabit, as they predict super- 
Eddington accretions. This is no longer true if one assumes 
that the Mbh - c c is independent of redshift, instead. Most 
of the other assumptions imply ~ 10 s - 10 10 Mq black holes, 
and accretion efficiencies of 0.01< e< 1. Our results suggest 
that at a given redshift, black hole mass is not strongly de- 
pendent on QSO bolometric luminosity, but a fixed value for 
the accretion efficiency is still a good fit to the data. 

These results are in agreement with those of 
iMcLure fc Dunlopl i|2(l(l 1 ; . In particular the latter measured 
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the masses and Eddington efficiencies of high-z black holes 
using data from the SDSS DR1, through modelling the QSO 
spectra. Their analysis, significantly different from that pre- 
sented here, results in Mbh and efficiency e values similar 
to those we obtained. Different relations between the black 
hole and dark halo masses differ almost by a scaling factor. 
Therefore, the trend observed in the Mbh — log(L) plot is 
the same irrespective of the halo density profile and Mbh - 
Mdmh;M B h - cr c relation. 
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